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Abstract. We have recently developed a geometric formulation of the stress field for an 
interacting quantum system within the local density approximation (LDA) of density 
functional theory (DFT). We obtain a stress field which is invariant with respect to 
choice of energy density. In this paper, we explicitly demonstrate this uniqueness by 
deriving the stress field for different energy densities. We also explain why particular 
energy densities give expressions for the stress field that are more tractable than others, 
thereby lending themselves more easily to first-principles calculations. 



Understanding a material's energetic response to strain is fundamentally cou- 
pled to understanding the physics of a broad range of phenomena, from surface 
reconstructions to piezoelectricity. For example, it has been demonstrated that 
knowledge of the spatial distribution of microscopic stress via first-principles cal- 
culations can help explain the onset of macroscopic polarization in a material [1]. 
Therefore it is important to understand the nature of stress at the microscopic 
level. 

Formally, the electronic contribution to the microscopic stress must be included 
via quantum mechanics. There have been numerous methods developed for comput- 
ing the quantum stress. (See Ref. [2] and references contained within.) The stress 
field, <r Q/ g(x), is a rank-two tensor field usually taken to be symmetric (torque-free). 
The divergence of a Q p(x) must equal the force field F a (x) of the system: 

F a = V p a aP . (1) 

(Note that the Einstein summation convention for repeated indices is used through- 
out the paper.) It is well-known that Eq. 1 does not uniquely define the stress field 
(regardless of whether a system is quantum-mechanical or classical) since one can 
add any tensor whose divergence is zero to 0"a/?(x) and still recover the same force 
field [3]. 

The volume-averaged or total stress T a p has been related to the energetic response 
of the system to uniform strain and to the integral over the stress field: 

Ta , = 9E_ 



= J a Q/3 (x)rf 3 x, 



(3) 



where e a p is a uniform scaling (strain) of the system. Nielsen and Martin showed 
that T aj 3 is a unique well-defined quantity which can be computed efficiently [2]. 
However Eq. 3 does not unambiguously define the stress field. One can write the 
total energy E as an integral over an energy density £(x). Therefore 

T ^ = [ d 3 x^-. (4) 

Comparing Eq. 4 and Eq. 3 could lead to a definition of the stress field as [4,5] 

^(x) = (5 ) 

However, the total energy can also be expressed many other ways: 

J d 3 xS + J d 3 xV a W a . (6) 

If variations of W a with respect to the dynamical variables (e.g. single-particle 
wavefunctions) vanish on the boundary of the system, then the minimum of E 
remains invariant with respect to choice of energy density. However, the energy 
density is changed by the presence of W , making the expression of Eq. 5 and Ref. 
4 dependent on the choice of energy density. Choosing a form of the energy density 
is referred to as specifying an energy density "gauge" , and hence the stress field of 
Eq. 5 is said to be not gauge invariant. 

We have recently developed a formulation of the quantum stress field which is 
invariant to choice of energy density, and have used it to successfully calculate the 
stress fields in periodic systems using LDA-DFT [6]. We define the stress field as 

de Q/3 (x) 

where e Q/ g(x) is the strain tensor field. This definition can be derived via a virtual 
work theorem [3] . We then use a relationship known in continuum mechanics that 
equates e a/ 3 with the Riemannian metric tensor <? a /?( x ) [7,8]. Therefore the stress 
field is 

2 5E 

where g and g alS are the determinant and inverse of g a /3, respectively. In order 
to calculate a a/3 , we first express the constrained energy functional for LDA-DFT 
[9,10] in a covariant manner via the principle of minimal coupling, and then com- 
pute the functional derivative using a procedure well-known in general relativity 



theory [6,11]. It should be mentioned that this prescription for obtaining a a p is 
not limited to LDA-DFT but can be generalized to derive the stress field for other 
DFT functionals as well as explicit many-body systems of quantum particles. 

We now demonstrate the invariance of our expression for the stress field with 
respect to choice of energy density. In DFT there are two terms in the total energy 
expression which are commonly subjected to energy density gauge transformations. 
The first is the single-particle kinetic energy which can be written as an integral 
over the so-called asymmetric kinetic energy density: 

= /-^E W<* (y/gs^drfi) d 3 x, (9) 

where the {0j} are the single-particle wavefunctions which obey the Kohn-Sham 
equations generalized for curvilinear coordinates [6]. Therefore we are consider- 
ing the total energy to be at a minimum. We consider here for brevity insulator 
occupation numbers. The variation of Eq. 9 can be written as 



SE k = sl\ /E Vgd^g^d^d 3 
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where we have used the covariant version of the divergence theorem in the second 
line. All variations in {0j} vanish on the boundary of the system, and we require 
5g a/3 = on the boundary. Therefore the surface term in Eq. 10 is zero. As a result 
the functional derivative required for evaluating Eq. 8 is 

5E^ 1 ^„ . d^g (1 



+ ^E/^w t P||. (id 

Note that this result is identical to what we would obtain if we wrote E k as an 
integral over the symmetric kinetic energy density: 

£ k = ^E/ y/gdafig^dfdHdPx, (12) 



and took the variation of this expression directly. Therefore the stress field is 
invariant with respect to choice of kinetic energy density. We mention that terms 
involving the variation of the wavefunctions with respect to metric vanish when the 



variation of the total ground-state energy is considered, since the wavefunctions 
obey the Kohn-Sham equations [6]. In other words, requiring 5E/5(p* = implies 



6g<#(y) dg<#(y) + ^ J 5<t>* 8<f#{y) *f J 5<pi 5g^(y) 



(13) 
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The Coulomb term describing the electrostatic electron-electron, electron-ion, 
and ion-ion interactions in the total energy functional is also commonly subjected 
to energy density gauge transformations. (This is entirely separate from transfor- 
mations of the U(l) electromagnetic gauge.) We can define -Ecouiomb as 

^coulomb = l - J ^/gpVd 3 x, (14) 

with 

p(x) = ^^(x-R l )-n(x), (15) 
i y9 

where Zj is the charge of the i-th ion located at position Rj, and n is the electronic 
charge density equal to J2i <f>i<f>i- The potential V can be computed from p via the 
Poisson equation: 

^pd a Lfgg af3 d p v) = -4np. (16) 
V 9 

The variation of -Ecouiomb is 

^coulomb = 5 {-i- J d a (V9~g a %v) Vd 3 x} 

= 5 {tl ^ d -V9 a %Vd 3 x^ -s{±Jd a (V^gg a %v) d 3 x} 
= 5 {tl ^ d « V 9 aP d(iVd 3 x} -8{±f Vy/ggpPdpVdS*} . (17) 

Since the variations of the potential and the metric vanish on the boundary, the 
surface term in Eq. 17 is zero. Therefore the functional derivative of -Ecouiomb with 
respect to metric is 

S -^coulomb _ J_ rv T 1 dy/g 

+ hJ^A* x > (18) 

where the electric field T a = —d a V. We would obtain this same result if we initially 
expressed -Ecouiomb as an integral over the Maxwell energy density: 




(19) 



This demonstrates that the stress field is invariant with respect to choice of elec- 
trostatic energy density. 

The non-local term involving the variation of the electric field with respect metric 
in Eq. 18 is unwieldy and is in general difficult to compute. To eliminate this term 
it is advantageous to choose an energy density that is the Lagrangian density in 
electromagnetism. Therefore we write -^coulomb as 



Now when -Ecouiomb is varied, the term containing 5V will vanish due to the Poisson 
equation. Variation of Eq. 20 also gives a term relating to the variation of p with 
respect to metric: 



As explained in the kinetic energy section, this variation (and those of other energy 
terms) will be multiplied by 5_E/<50*(x). The Kohn-Sham equation insures that 
variations of E with respect to 0* vanish, and the ionic charges are fixed, thereby 
removing the need to evaluate Eq. 21. 

In conclusion, we have demonstrated explicitly that our formulation for the quan- 
tum stress field within DFT is invariant with respect to choice of energy density. 
Therefore the stress field is a well-defined object that can be computed via first- 
principles to help gain a microscopic understanding of stress-mediated phenomena 
in complex materials. 

The authors wish to acknowledge R. M. Martin and D. H. Vanderbilt for helpful 
discussions. This work was supported by the Office of Naval Research under grant 
number N-00014-00-1-0372 and the Air Force Office of Scientific Research, Air Force 
Materiel Command, USAF, under grant number F49620-00- 1-0170. 



1. N.J. Ramer, E.J. Mele, and A.M. Rappe, Ferroelectrics 206-207, 31 (1998). 

2. O.H. Nielsen, and R.M. Martin, Phys. Rev. B. 32, 3780 (1985). 

3. L.D. Landau and E.M. Lifshitz, Theory of Elasticity, 3rd ed. (Butterworth- 
Heinemann, Oxford, 1986), pp. 1-7. 

4. A. Filippetti and V. Fiorentini, Phys. Rev. B. 61, 8433 (2000). 

5. N. Chetty and R.M. Martin, Phys. Rev. B. 45, 6074 (1992). 

6. C.L. Rogers and A.M. Rappe, Phys. Rev. Lett, (submitted 2000) cond-mat/0006274. 

7. Y.C. Fung, Foundations of Solid Mechanics, (Prentice-Hall, Englewood Cliffs, 1965), 




(20) 




(21) 



REFERENCES 



pp. 90-92. 



8. L. Mistura, J. Chem. Phys. 83, 3633 (1985); Inter. J. Thermophys. 8, 397 (1987). 

9. P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

10. W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

11. L.D. Landau and E.M. Lifshitz, The Classical Theory of Fields, 4th ed. (Pergamon, 
Oxford, 1975), pp. 270-273. 



